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ABSTRACT 

One of the persistent complications in searches for transiting exoplanets is the 
low percentage of the detected candidates that ultimately prove to be planets, 
which significantly increases the load on the telescopes used for the follow-up 
observations to confirm or reject candidates. Several attempts have been made 
at creating techniques that can pare down candidate lists without the need of 
additional observations. Some of these techniques involve a detailed analysis 
of light curve characteristics; others estimate the stellar density or some proxy 
thereof. In this paper, we extend upon this second approach, exploring the 
use of independently-calculated stellar densities to identify the most promising 
transiting exoplanet candidates. We use a set of CoRoT candidates and the 
set of known transiting exoplanets to examine the potential of this approach. 
In particular, we note the possibilities inherent in the high-precision photometry 
from space missions, which can detect stellar asteroseismic pulsations from which 
accurate stellar densities can be extracted without additional observations. 

Subject headings: Stars: Planetary Systems, Techniques: Photometric, Stars: 
Fundamental Parameters (density) 
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Introduction 



One of the principle goals of exoplanet searches using transits is the discovery of 
temperate terrestrial exoplanets. Their detection presents a tremendous challenge - 
particularly for the radial velocity observations that are typically used to confirm a 
transiting exoplanet candidate. Even the confirmation of transiting Hot Jupiters can 
confound observers, consuming and ultimately wasting telescope reso urces. Transit surver ys 
typicall y identify 10-20 cand idat es for each pl anet (STAR. 



OGLE flUdalski et al. 



2002a 



20031 ). HAT flBakos et al 



, WASP fiPoUacco et al. 



20061), 



2004 )). although, in the case of 



CoRoT, t he ratio of top 



one-half ( ICabrera et al. 



prior ity candidates from that mission is closer to one-third or even 



20091 ). This demonstrates the success of a decade-long effort to 



develop tec hniques to elminate tran sit candidates with further observa tions. Out-of-ecli pse 



variations ( ISirko fc Pacyhski 



stellar den sity diagnostics (ISeager fc Mallen-Ornelas 



2003h. V-shaped vs. U 



diagnostic fiTingley fc Sackett 



shape d transits ( jPont et al. 



20051 ) 



20031 ). and the Tingley-Sackett 



20051 ). to name a few, have all contributed. Combined, these 



technique have enabled CoRoT's relatively high success rate compared to, for example, 
the first OGLE run, which found 2 planets among 59 equally-prioritized candidates 



(Udalski et al. 



( jLeger et al. 



2002al Jbr). Even so, many candidates defy easy characterization. GoRoT-7b 



20091 ). to use the lightest known__tran^itiiig 



more than 100 observations with HARPS ( lOueloz et al. 



exop 



anet as an example, required 



20091 ). totalling over 70 hours of 



observing time. While the host star is fairly active in this case, thus complicating attempts 
to detect the radial velocity signal of the planet, the signal of GoRoT-7b was more than 
an order of magnitude larger than what one would expect for an Earth twin (3.3 m/s vs. 
0.1 m/s). With the anticipated influx of candidates consistent with terrestrial exoplanets 
from both CoRoT and Kepler (and potentially future missions, such as PLATO), as well 
as from ground-based searches for transits across low-mass stars, any technique that could 
eliminate even a few terrestrial-class candidates would potentially save an enormous amount 
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of valuable telescope resources. 



One avenue of candidate screening that has not be e n ful 



diagnostic first presented by 



Seager fc Mallen-Ornelas 



y explored is the density 



(120031 ). Here, the authors 



demonstrated that it was possible to extract stellar densities by fitting the transit light 
curve with a simple trapezoid, which we will refer to in this paper as psMO- At the time 
that paper was written, the required photometric measurements (5 mmag precision with 
5 minute sampling covering two transits) was not at all typical for transit candidates - in 
fact, it was quite unusual. This paper was the first in what has turn ed into a long discus sion 



on wh at exactly could be done with detailed analyses of transits. 



Tinglev fc Sackett 



(120051 ) proposed a technique to identify the best transit candidates in a sample using only 



transit periods, durations, and depths and inc luded a discussion of the impact of unknown 



eccentricity on their technique. 



Barnes 



(120071) and 



Burke 



(I2OO8I ) later expanded upon this 



secondary analysis to conclude that transits are more likely to occu r in the elliptic a. 



and this results in highe r yields from transit surverys, res pectively 
revisited one part of the 



Sozzetti et al. 



orbits 



Seager fc Mallen-Ornelas 



(I2OO7I ) 



possible to obtain stellar densities from transit fits and that t 
the derivation of stellar parameters from evolutionary tracks. 



(120031 ) derivation, stating that it was 
lese could be used to aid in 



Ford. Ouinn fc Veras 



(120081) 



showed that it was poss i ble to estimate unknown eccentricites through an analysis of the 



transit shape; 



Kipping 



(I2OO8I ) did something similar, d eriving an enti r ely ne w equation 



Yee &: Gaudi 



Kipping 



(l2008h w orked to 



that involved no potentially erroneous approximations, 
estimate period for transit candidates with but a single transit, 
the accuracy of different equations for t ransit duration, assessing the effect of incorrectly 



( 120101 ) compared 



assuming circular orbits. 



Brown 



( I2OIOI ) discussed the use of stellar densities from transits 



to restrict stellar radii, which are crucial in determining planetary radii. 



In this paper, we combine and expand upon these ideas and propose a slightly different 
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technique to identify good transit candidates: comparing densities calculated from transits 
to densities calculated by other, independent means, such as J — K colors, spectra in 
combination with evolutionary tracks, and asteroseismology. It is our feeling that this 
will be of particular interest for shallow candidates discovered using ultra-high precision 
photometry from space-based mission such as CoRoT, Kepler and potential future missions 
such as PLATO for which the expected radial velocity signal would be weak and difficult 
to detect. To evaluate this technique, we use both the set of CoRoT candidates and the 
known transiting exoplanets to demonstrate that this diagnostic has the potential to be 
very effective in the pre-selection of the best candidates. To aid in this process, we derive a 
new equation for calculating stellar densities using other, more tangible transit parameters 
often available in papers that quote neither a/i?* nor the stellar density. This equation 
includes orbital parameters (eccentricity and angle of periastron), which make it possible to 
analyze the impact of an unknown orbita l ecce ntricity on the measured s tellar density. It is 



based on the equations in 



Sackett 



fll999h and 



Tingley &: Sackett 



(I2OO5I ) and may also be 



used as a consistency check of fitted transit parameters. 

In § [21 we derive the equations that will be used for the analysis. In § [3l we discuss the 
independent density measures that compliment the stellar densities from transit parameters. 
In § m we apply these equations to those transiting exoplanets for which the necessary 
parameters have been published. In § |5l we examine the use of the psMO combined with 
the densities derived from J — K colors on the CoRoT candidates. In § |6l we analyze the 
impact of eccentricity on this technique. Lastly, we discuss our conclusions in § [71 



2. Equations 



Calculating the stellar densities from transit parameters is not difficult. The most 
straightforward way begins with the equation for the density of a spherical star: 
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3M, 



(1) 



where is the density of the star and is the mass of the star. Then, one can use Kepler's 
3rd law to substitute out M^, needing only the rather safe assumption that ^ Mp, 
obtaining an equation for the density of the star based on the transit parameter a/i?^ {Pti}'- 



Pti 



where G is the gravitational constant and T is the orbital period of the planet 



(2) 



(ISeager fc Mallen-Ornelas 



2003 



Sozzetti et al. 



20071 ). Despite the absense of parameters 



associated with orbital eccentricity, this equation does yield accurate stellar densities from 
transits of planets in elliptical orbits. Either the eccentricity and argument of the periastron 
are supplied from the analysis of the radial velocity observations and used during the transit 
fit or the radial velocity observations and the photometry are fit simultaneously. 

It must be pointed out that pti is distinctly different from psMO, despite the fact that 
each are derived in the same article, pn uses parameters extracted from a full transit fit with 
limb darkening to estimate the stellar density, while psMO uses a trapezoid fit to estimate 
the stellar density, neglecting limb darkening. Moreover, psMO assumes circular orbits, while 
Pti in principle takes into account eccentricity and the argument of the periastron. 

It is possible to use a/i?^ as a free parameter when fitting a transit, but many papers 
describing transit observations do not include it. Moreover, one cannot intuit whether the 
value makes sense from looking at a transit, which is clearly a problem in a few cases. 
Therefore, it is instructive and useful to derive equations for the stellar density and for 
a/i?^ that depend on transit parameters t hat a re tangible and more commonly published. 



Starting with Eq. 7 of 



Tingley fc Sackett 



(itself derived from equations in 



Sackett 



(|l999|)): 



20051) for the duration {th) of a planetary transit 
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^ , v^r^e2 / a2(l-e2)2 cos2 i ( T \ - 

where i?p is the radius of the planet, e is the planet's orbital eccentricity, uj is the argument 
of periastron, and i is the orbital inclination. It assumes that the mass of the planet is 
negligible, the star-planet separation during transit is much greater than the stellar radius 
(which is not always the case), and the orbital velocity is constant during the transit. From 

— 1/3 

inspection, it is evident that the right-hand side of above equation goes as RM* ' , which 
is proportional to p^^^^. If we solve for M}^^ /R^, cube everything and multiply by 47r/3, we 
obtain another equation for the density {pt2) of the parent star, based on more fundamental 
parameters than a/R^: 



^*^ = ^l l-e2 ) 



where the ratio of the radii k = Rp/R^,, the impact parameter b = Qa cos i/ R^, and 

1 - 

Q = 7—^' (5) 
1 + e sm a; 

which equals 1 for circular (e = 0) orbits. One can also solve the same equation for a/i?^, 
getting: 



f = ^^Vil + ky-b'- (6) 

-K* 7rri4vl — 

Plugging this equation into Eq. |2] yields the same results as given in Eq. H] . As this 
equation is invertible, transit durations can be calculated from a/i?^ as well. This makes it 
possible to check the quahty of a fit in cases where all of the pertinent parameters (T,ri4, 
a/i?^, b, k, e, u) are given. Additionally and significantly, this equation makes it possible 
to evaluate the potential impact of an unknown orbital eccentricity on densities determined 
from transit parameters (referred to generically as transit densities or pt hereafter). 
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With the se two equation s in hand, we can assess if the assumption of a ^ i?^ is an 



issue - indeed 



Kipping 



( 120101 ) specifically calls into question the accuracy of this equation - 



but given the level of uncertainty in transit parameters, it ultim ately make little dis cernable 



2009aJ). Even 



difference. HAT-P-7b has one of the smallest a/-R^s, only 3.82 ( IWinn et al. 
for this case, the impact on the density is very small: pti = 0.217 and pt2 = 0.204 - only 6%, 
which is approximately one-third the size of the measurement errors and not an atypical 
difference between pa and pt2 for the transiting exoplanets for which both values are 
available. However, planets may someday be discovered transiting very closely to evolved 



or very hot stars - for such cases, this assumption might very well prove problematic. 



3. Independent Stellar Density Measures 

In order to accomplish our stated goal of evaluating transit candidates by their stellar 
densities, it is necessary to have an independent measure of stellar density to compliment 
the transit density. In the following analyses, we pursue three different possibilities: 1) from 
published values of stellar radii and masses not based on transits, 2) from J — K colors, and 
3) from asteroseismology. As a group, these will hereafter be referred to as stellar densities. 



3.1. Stellar Densities from Published Stellar Masses and Radii 



It is a trivial matter to calculate stellar densities from the published stellar masses 



and radii (referred to hereafter as pspec)- These are gene rally derived 



Sozzetti et al. 



rom spectral 



(120071 ) for a 



analysis combined with modeled evolutionary tracks (see 
detailed discussion). Spectral analysis yields a variety of information, such as temperature, 
metallacity, logg (where g = GM^,/Rl), which can then be used with either a true distance 



measure (generally from Hipparcos) or evolutionary tracks to obtain the stellar radii and 
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(120071 ). it is possible to 



densities, using p*/ Pq = q^/qqE^. According to ISozzetti et al 
use logg to constrain the stellar parameters from the evolutionary tracks, but surface 
gravities from spectra generally are not very precise, leading to highly uncertain stellar 
radii. They then proposed using transit densities to improve the determination of stellar 
parameters. This would suggest that at least some groups have used the transit densities 
to help deterrnine th eir stellar parameters; while it was not done before the TrES-2 paper 



( ISozzetti et al 



20071 ). it has in the intervening years become essentially universal. Therefore, 



we caution that for some published transiting exoplanets, this density {pspecj may not be 
truly independent of the transit density. 



3.2. Stellar Densities from J — K colors 



Allen's Astrophysical Quantities (jCox 1 120001 ) lists J — K colors and stellar masses 
and radii as a function of spectral type for main sequence stars. While these values are 
very general, being independent of factors such as age and metallicity, they may offer 
a reasonable approximation of stellar density, at least for the purpose of identifying 
less-than-ideal candidates. The advantage of this m easure is that J — K colors are readily 



20061 ). which means 



available for most stars from the 2MASS database (iSkrutskie et al. 
that no additional observations or analysis are necessary - simply look up the J ~ K color, 
then interpolate among the values from Allen's Astrophysical Quantities to arrive at the 
density (pjk) considering only main sequence stars. This is illustrated in Fig. [H which 
plots the stellar density versus the J — K color. 

Apparently, this measure has several weaknesses. First, it assumes that all stars are 
main sequence stars. Obviously, this is not always the case; stars evolve, and once off the 
main sequence the resulting Jk densities would be incorrect. Second, while J — i^' is a very 
red color, stars residing in dust-rich environs can experience enough interstellar reddening 
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to have a significant effec t on their resu 



excess E{J — K) = 0.24 (IBonomo et ah 



ting 



density. For example, CoRoT-lOb has color 



2OIOI ). which leads to a difference in log(p^/p0) of 



about 0.5 dex - or a factor 1.6. 



Technically speaking, the 2MASS J and K colors are slightly different than the 



Bessel and Brett 



Bessel fc Brett 



(119881 ) system used in Allen's Astrophysical Quantities. 
Therefore, to be completely rigorous, a conversion from the one system to the other ought 
to be performed. However, such a conversion results in a very small change in stellar density 
- only a few hundreds of a solar density. Considering the magnitude of the generalizations 
involved in this density measure, this difference can be neglected. 



3.3. Asteroseismic Densities 



Another technique exists that can reliably obtain stellar densities, though we cannot 
exploit it yet in this work: asteroseismology. By measuring the frequencies of the acoustic 
modes and performing some model-fitting, it is possible to get a good measure of the 
stellar density from the ay e rage frequency spacing between these modes in the asymptotic 



region ( iVandakurov 



1967 



Tassoul 



19801). It is only in the last few years that these 



measurements hav e become possible, generally through the radial yelocity fingerprint of 



these oscillations (le.g. Kjelsen. Bedding, fc Christensen-Dalsgaard 



(I2OO8I)). Recently, 



however, circumstances have changed: the Kepler satellite mission has the capacity to 
obtain photometric time series with the necessary precision to reveal stellar oscillations, as 
would proposed missions in the future such as PLATO. This is potentially very interesting: 
the photometric precision necessary to discover transiting terrestrial exoplanets should also 
be sufficient to detect oscillations, thereby offering a sensitive density measure without 
additional observation. 
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Fig. 1. — J — K colors and spectral types vs. stellar density. This figure shows how the 
stellar density varies with J — K color and spectral type for the rnain s equence. Diamonds 



show the actual values from Allen's Astrophysical Quantities (ICox 



20001), while the line and 



asterisks show the interpolated values (a fourth-degree polynomial fit) used to calculate the 
J — K stellar densities used in this paper. 
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4. Application to Known Transiting Exoplanets 



We performed an extensive literature search to gather the transit parameters for as 
many transiting exoplanets as possible. Table 1 lists the planets and their periods, along 
with their corresponding pn and pt2 (calculated from the available transit parameters) 
and Pjk, when possible - some faint (or very crowded) candidates do not appear in the 
2MASS databas e. Asteroseismic densities are much less common, however: only a single 



recent paper by 



Christensen-Dalsgaard et al. 



(j201o| pubhshes asteroseismic densities (pa) 



and then only for the three stars known to harbor exoplanets in the Kepler field before the 
mission began: HAT-P-7, HAT-P-11, and TrES-2 - the last two with preliminary stellar 
parameters only. These are also included in Table 1. 

Fig. [2] shows the ratio of the transit densities and the stellar densities from masses 
and radii, with pti on the top and pt2 on the bottom. These figures demonstrate that, 
for nearly all of the candidates, the two values agree within the errors bars - even the 
eccentric ones. If anything, they agree far more than we would expect from the size of the 
error bars - evidence that the practice of using transit densities to constrain model fits of 
stellar parameters is widespread. Only four planets seem to have deviations greater than 
la from a density ratio of one for either pn/pspec or pt2/Pspec- HAT-P-7b, CoRoT-llb, 
CoRoT-7b and 0GLE-TR-211b. The first case is resolved w hen the transit parameters 



are taken fro m the same referenc e as the stellar parameters (jPal et al. 



a later paper (I Winn et al. 



2008|), instead of 



2009af ). which contains better photometry and consequently 
more precise transit parameters but does not repeat the spectral analysis of the parent 
star. This presents a clear case where the derived transit density influences the stellar 
parameters. The second offers no easy explanation and requires further study. The last two 
do not have sufficiently precise photometry for confident parameter estimation: CoRoT-7b 
has an extremely shallow transit and it has been further hypothesized that magnetic 
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activity or transit timing variations may affect the shape of the transit f lLeger et al.ll2009l ) 
and 0GLE-TR-211b has not been observed well enough - several different transits are 
spliced together to create the final light curve, introducing undefined uncertainties. This 
emphasizes the fact that good photometry without significant unknown systematics is 
essential for the use of this technique. 

As can be seen in Fig. |3l the details of the distribution of the density ratios is 
interesting. The region is largely triangular, with a narrower distribution for deep transits 
and wider for shallower ones. This apparent dependency could arise from two causes: 
underestimations of pu and pt2 or overestimation of Pjk- It is unlikely that bad transit 
parameters are the source of this problem: while shallower transits would tend to have less 
precise transit parameters, the precision of the transit parameters for the CoRoT planets is 
more limited by uncertainties in the limb darkening parameters than photometric precision 
and they exhibit the same tendency. Therefore, we suspect that pjx is the source. Most of 
the planets - but of course not all - discovered are approximately Jupiter-sized. Therefore, 
some fraction of shallower transits correspond to larger (evolved) stars. The Pjk^ derived 
are based on the assumption that the stars are still on the main sequence. If the J — K 
colors accurately represent the densities of main sequence stars, some of the larger stars will 
be more evolved stars, rather than higher mass, explaining the presence of unexpectedly 
low density ratios, leading to the observed increase in the spread in density ratios as transit 
depth decreases. The numbers actually bear this out; the four exoplanets with the lowest 
density ratios (HAT-P-4b, HAT-P-13b, HD 149026b, and HD 17156b) have orbit stars with 
effective temperatures between 5650 and 6150 Kelvin, but R ~ 1.5Rq and logg ~ 4.2, 
when ~ 4.5 is more typical for main sequence stars with these temperatures - indeed, HD 
149026b is classified ClS db subgiant. 

Asteroseismic densities may be more accurate than Pjk, the few known examples 
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exhibiting little scatter around the expected density ratio of unity. While the current 
sample is too small for conclusions, we note that additional observations are not necessary 
to obtain the asteroseismic density, only some analysis. Therefore, a comparison 
between asteroseismic densities and transit densities, when possible, should prove to be 
an extremely valuable diagnostic tool in evaluating transit candidates. This may be 
particularly useful for candidates for terrestrial exoplanets in multi-planet systems - some 
groups have hypothesized that terrestrial planets m ay favor low-eccentric ity orbits to 



maintain stability, both with g iant planets present ( iPilat-Lohinger 



( iMann. Gaidos. &: Gaudi 



20091 ) and without 



20101) ■ This topic is still being debated, however, and does not 



address single planet systems al all. 



5. Application to CoRoT Candidates 



The set of CoRoT candidates offers an excellent opportunity to test this technique. 
Many of these can didates are listed in p apers describ ing the results of ind ividual CoRoT 



runs - e.g. IRa02 (iCarpano et al. 



20091 ) and LRcOl ( jCabrera et all 120091 ) - but others 



are not yet published. Ideally, we would want to use transit densities calculated from 
detailed transit fits and compare them to densities from spectra, but this information is 
not currently available for more than a handfuU of candidates. However, as part of the 
array of information used for ranking candidates, the CoRoT Science Team determines 
both psMO (derived from a trapezoid fitting that neglects limb darkening) and pjx- Each 
of these techniques is somewhat less than ideal, but without spectral information for each 
candidate, it is impossible to use detailed transit fitting including limb darkening and stellar 
evolutionary tracks to get better densities. They are apparently sufficient, however, to test 
the potential of density for candidate evaluation. As can be seen in Fig. |U these simple 
density measures work remarkably, perhaps even surprisingly, well. While not all of the 
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Fig. 2. — The left-hand figure shows the density ratios for the known transiting exoplanets 
for which the necessary information is available. Density ratios are created using spectra 
(p^) and either pt\ (top) or (bottom). Notice how consistently these values are centered 
around 1, particularly for pt\ - clear evidence of the widespread use these densities have 
in the final determination of stellar parameters. The right-hand figure reiterates this with 
a comparison of the density ratios calculated using pti and pt2 for the subset of transiting 
exoplanets for which all the necessary parameters are published. Note that the values for pt\ 
and pt2 are very similar in most cases, even for planets with high eccentricities, verifying that 
Pt2 is accurate. Those that appear to be outliers actually aren't - the errors bars are large 
enough to include the anticipated result. Those that deviate from ratios of 1 but remain 
close to the overplotted fine probably have problems obtaining the stellar parameters from 
the spectra (and models), while those that are displaced upwards and downwards typically 
have problems with the photometry. 
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Fig. 3. — Ptransit/ P-k VS. transit depth for known exoplanets. This figure shows the density 
ratios (in this case, Ptxj Pjki when available, otherwise Ptij Pjk) for all the known transiting 
exoplanets (squares) and for the three cases for which the asteroseismic density {pti/ Pa or 
Pnl Pa is known (diamonds)). The depth is approximated as the square of the planet-star 
radius ratio. The known exoplanets have density ratios distributed in a way that appears 
to be asymmetric relative to the ideal Ptransul P-k = 1 value, likely due to stellar evolution 
(shaded region). The limited sample of asteroseismic densities appears to be more reliable, 
agreeing well with the transit densities regardless of transit depth. 
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Fig. 4. — Ptransit/ P* VS. transit depth plot for the CoRoT candidate hst. Diamonds are the 
confirmed exoplanets, while large circles are the highest priority CoRoT candidates and small 
circles lower-priority CoRoT candidates. For these, psMO is the transit density estimated by 
the trapezoid method and pjk is the stellar density. The shaded region is the same as in 
Fig. m Notice how the CoRoT planets fall in the same shaded region. Also notice that the 
vast majority of candidates have very low density ratios, while very few have ratios higher 
than 1. This is because candidates in both giant stars a nd significantly blended stars will 



have low density ratios (ISeager fc Mallen-Ornelas 



20031 ). while only eccentricity or errors 



can increase density ratios. This demonstrates how effectively density ratios can be for 
screening transit candidates. It also shows that psMO is a reasonable proxy for pn and pt2- 
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CO (degrees) 



Fig. 5. — Density Ratio vs. eccentricity and angle of periastron. This figure shows how 
the density ratio varies as a function of angle of periastron and eccentricity. The lines from 

outermost to innermost depict eccentricities going from 0.9 to in steps of 0.1. Notice how 
the angle of periastron produces density ratios greater than 1 over a larger range of values 
than it does for ratios less than 1. 
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0.01 0.10 1.00 10.00 100.00 
Pt2(e,cj)//),2(e = 0) 

Fig. 6. — Density ratio vs. probability and eccentricity. This figure shows the probabihty of 
measuring a particular density ratio as a function of eccentricity. The lines from innermost 
to outermost depict eccentricities going from 0.9 to 0.1 in steps of 0.1. Notice how the 
probability that the density ratio is greater than 1 increases with eccentricity. 
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Fig. 7. — Density ratio probabilities vs. period. This figure shows the probabihty that a 
transiting planet with unknown e and uj will produce a given density ratio for three different 
period regimes, using the observed distribution of exoplanet eccentricities. The dashed line 
is for planets with periods less than 5 days, while the dotted hne is for planets with periods of 
5 to 100 days and the solid line is for planets with periods longer than 100 days. Notice how 
the peak at a density ratio equal to one decreases as period increases, while the chance for 
extreme density ratios increases. Please note that the lines have been smoothed for clarity. 
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Fig. 8. — Density ratio probabilities vs. eccentricity. This figure shows the probabihty that 
a given eccentricity with an unknown periastron angle will have a transit density lower than 
the listed value. The eccentricity ranges from 0.1 to 0.9 in steps of 0.1. From this plot, it 
is clear that high eccentricities strongly favor transit densities greater than 1, while low-to- 
moderate eccentricities have a significant probability of producing comparatively low density 
ratios. Arbitrarily choosing a density ratio cut-off of > 0.5 would falsely de-prioritize ~ 11% 
of transits with an eccentricity of 0.7, , ~16% of transits with an eccentricity of 0.4 and no 
transits with an eccentricity less than about 0.22, while including relatively few astrophysical 
false positives. It should also be mentioned that low eccentricities < 0.3 are more common, 
particularly for exoplanets with periods less than 5 days - which transit searches are most 
likely to find. 
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planets line up along the line describing the density ratio psuo/p* = 1, niost of the planets 
do however fall in or are very close to the same well-defined region as the planets shown in 
Fig. 131 We note that the follow-up process is incomplete and that some of the candidates in 
that region may in fact be unconfirmed planets. 

It is also clear in this figure that the majority of candidates have very low density 
ratios. This can be primarily attributed to two causes beyond the limitations of psMO- 
giant stars and blends, which are eclipsing binaries whose light is accompanied by a bright 
third star, causing de ep eclipses to be diluted and thu s appear like planetary transit. As 



already mentioned in 



Seager &: Mallen-Ornelas 



( 120031 ) ■ both blends and giant stars will 



have unusually low values of psMO ^ the giants stars because their density is truly much 
lower than main sequence stars and the blends because the inclusion of the light of the third 
star leads the trapezoid fit to converge on a solution that is larger than the bright third star 
and with a higher impact parameter. This overestimation of the radius leads in turn to an 
underestimation of the stellar density of the third star of up to 50%. 

The known planets in this figure support the hypothesis that the distribution of 
densities ratios is due more to the insufficiency of pjx than the transit densities. The 
precision of CoRoT data is extremely high, so the errors in the transit parameters are 
dominated by uncertainties in the limb darkening rather than limitations of the photometry. 
Indeed, despite the fact that psMO is built upon more risky assumptions than pn and pt2 
and therefore presumably less reliable, the distribution of the psMO densities for the CoRoT 
planets are essentially indistinguiable from the distribution of the pn and pt2 densities for 
all exoplanets (including the CoRoT planets) - while the individual values are different, 
the distribution is not. The only exception is CoRoT-7b, for which the extracted transit 
parameters are likely distorted by magnetic activity or transit timing variations. This 
demonstrates that psMO is actually quite useful in and of itself (particularly for more central 
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transits) given the errors inherent in the extracted transit parameters, which may or may 
not include impossible-to-define uncertainties in hmb darkening coefficients. 

6. Impact of Eccentricity on Densities from Transits 

As the transit density is a function of the orbital eccentricity and angle of periastron, 
the danger in using this diagnostic blindly is the high probability of a significantly non-zero 
eccentricity - most planets have significant eccentricity, particularly among planets with 
periods greater than 5 days or so. To evaluate the value of this technique properly, it is 
necessary to evaluate the infiuence of non-zero eccentricities. 

It is a relatively straightforward matter to use the equation for pt2 (which directly 
includes eccentricity while the equation for pn does not) to explore this dependency. By 
comparing the response of pt2 as e and uj are changed to pt2 with e = 0, we can observe the 
effect of unknown eccentricity on the Pt2(e, a;)/pt2(e = 0) density ratios. This is shown in 
Fig. [5l These can readily be converted into the probability of observing a particular density 
ratio as a function of e, if u is unknown. This is shown in Fig. [61 Here, we see that, as 
eccentricity increases, the probability of observing a density ratio greater than 1 increases. 
This is the case because the probability that a transit occurs goes as the inverse of the 
star-planet separation, which therefore more likely that a transiting planet will transit 
when it is close to its parent star, rather than far away. The closer the planet is to its 
parent star, the higher its orbital velocity relative to its hypothetical e = velocity. Higher 
velocities lead to shorter transit durations and thus higher transit densities. The higher the 
eccentricity, the more pronounced this effect: 93.9% of all orbits with e = 0.9 will produce 
a density ratio greater than 1. This can then be combined with the observed distribution 
of planet eccentricities as a function of period (Fig. [7j), to give an idea of what fraction 
of planets with a given density ratio and period are planets. It is important to remember 
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that this sample is comprised almost entirely of giant planets, however - the distribution of 
eccentricity as a function of period might be different for terrestrial planets. 

One final consideration must be remarked upon: while density ratios larger than 1 
are more common for planets of unknown eccentricities, very low density ratios will also 
occur, albeit rarely. These will exhibit very distinctive light curves that display very long, 
fiat-bottomed transits. A density ration of 0.01 would result in an increase in the duration 
by a factor of ^1/0.1 = 4.6 over that expected for a circular orbit with the same period. 
Given that planets in circular orbits with periods of a year would have transit durations 
over 12 hours, such planets could have transit durations up to several days, but would stand 
out among false positives with similar density ratios. To the best of our knowledge, only 
eclipses of giant stars could produce similar events and these can be quickly identified with 
spectra due to their low surface gravity. For the interpretation of transit candidates, this 
means that cases with low pt/ Pjk (or Pt/ Pa-, if available) should first be studied to see if 
the prospective parent star is a giant. If not, the probability that the candidate is a planet 
remains very low (see Fig. [8]) and should thus be assigned a very low priority. 

To apply this method, we recommend selecting a lower limit for the density ratio as 
a function of the period: only planets with periods > lOOrf are currently known to have 
a chance to have very extreme eccentricities; in fact, the maximim po ssible eccent r icity 



Deeg et al. 



mm)- 



dereases with period: e < 0.35 log(T/l(i) + 0.24 (based on figure 4 from 
Accordingly, no known 'hot' planets (T < 5c/) have an eccentricity higher than 0.265. With 
this in mind, we can therefore conclude, per Fig. [HI that all candidates with density ratios 
lower than 0.45 can be strongly de-prioritized. Similar 'secure' cut-offs could be established 
for candidates with longer periods, which could potentially possess larger eccentricities. 
Additionally, 'semi-secure' cut-offs could be implemented, which would have small residual 
probabilities of falsely de-prioritizing very eccentric planets. However, should pjk be used 
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to create the density ratio, it is perhaps better to be more circumspect given this density 
measure's difficulties with evolved stars. 

7. Conclusions 

Ratios of densities derived from transit parameters over those from J — K colors 
(or some other independent measure, such as asteroseismology) are a useful tool for the 
identification of transit exoplanets out of a much larger set of candidates. Even rudimentary 
density measures (trapezoidal and J — K densities) appear to identify exoplanets effectively 
among the CoRoT candidates. More precise density measures from detailed transit fitting 
and stellar spectra will be even more effective; these, however, require additional spectral 
observations. Given data of sufficient quality, densities from asteroseismology might be even 
more effective, although at this time the sample of systems with transiting exoplanets and 
asteroseismic densities is too small to draw a definitive conclusion. Significant eccentricities 
will impact the measured density ratio in the absense of any knowledge of these orbital 
parameters; however, they will tend to increase the density ratios of discovered transiting 
exoplanets, while most false positives have density ratios lower than typical planets. It is 
important to remember that any tool can be misused; the one describe in this paper is 
no exception, but it does represent an improvement on other techniques currently used to 
prioritize transit candidates using transit parameters. Even a 'secure' density cut-offs as 
described in the previous section may de-prioritize some true exoplanets due to any number 
of causes: extreme planet or star characteristics or otherwise inaccurately measured density 
ratios, for example. However, when used in conjunction with other techniques (ellipsoidal 
variations, presence of secondaries, etc.), density ratios can help identify the candidates 
which are most likely to be exoplanets and which are most likely to be false positives, 
increasing the efficiency of efforts to confirm or reject planet candidates. 
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